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Abstract. Unitarity corrections to several electromagnetic processes in collisions of relativistic heavy nuclei 
are considered. They are due to the unitarity requirement for the S-matrix and correspond to the exchange 
of light-by-light scattering block between colliding nuclei. We obtain improved results for the corrections 
to e + e~ and p + p~ P a i r production, as well as new results for unitarity corrections to the production of 
photons via virtual Compton and virtual Delbriick scattering. These corrections can be numerically large; 
e.g., the p + p~ pair production cross section is reduced by about 50 % and the nuclear bremsstrahlung by 
about 15 -r 20 %. 



1 Introduction 

The subject of this paper are so-called unitarity correc- 
tions, which form a conceptually interesting class of cor- 
rections for quantum electrodynamic (QED) processes. 
The unitarity corrections come from the unitarity require- 
ment for the S'-matrix and are relevant for processes in 
which the lowest order amplitude is large, i.e., in which 
a large number of photons and/or real electron-positron 
pairs are typically produced. One class of such processes 
is bremsstrahlung and lepton pair production in ultra- 
relativistic heavy-ion collisions. 

In order to put the current investigation into perspec- 
tive, let us briefly recall that heavy-ion collisions definitely 
concern matter under extreme conditions. The impact pa- 
rameters p relevant for all cross sections under study in 
the current paper are in the range of the electron Compton 
wavelength, 

p ~ A e = — = 386fm. (1) 
mc 

Now, the typical electric fields, as seen in the laboratory 
frame, generated by a particle moving at a speed charac- 
terized by a relativistic factor 7, are of the order of 



£( Cl) ~ TT 7 



Z e 7 m 2 c 2 
7 = £ cr Za 7 . 



he he 



(2) 



where £ a is Schwinger's critical field strength. The quan- 
tity Za^i assumes values in the range of ~ 60 for RHIC 



and ~ 1800 for the LHC, illustrating that the flash field 
accompanying the nuclei may well exceed the critical field 
by several orders of magnitude, for very small spatial re- 
gions and very small times. 

In view of the huge pair production and photoproduc- 
tion cross sections encountered under these conditions, it 
is natural that the importance of unitarity corrections has 
been recognized for the first time within the production 
of electron-positron pairs in collisions of heavy nuclei (see 
Rcf. pQ and the reviews |2 , 3J ) . The unitarity correction 
for the one electron-positron pair production process 



Z\ + Zi — ► Z\ 



e ' e 



(3) 



has been calculated in Ref. [I] and found to be about 
3 -7- 4 %. Estimates of unitarity corrections for the p + p~ 
single-pair production process 



Zi + Z 2 — > Zi + Z 2 + p n~ 



(4) 
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have been obtained in Ref. [5]; in this case the unitarity 
correction is found to be large (~ 50 %). 

Let us consider these corrections conceptually using 
the process (]3]) as an example. In this case the lowest or- 
der Feynman diagram is represented by Fig. [TJ while the 
diagrams of the type depicted in Fig. [5] correspond to the 
unitarity correction. These diagrams include blocks of vir- 
tual light-by-light scattering via an electron loop, whose 
imaginary part corresponds to the production of electron- 
positron pairs by the Cutkosky rules. 

For Z\ Z2 a 3> 1 and 7 ^ 1, it is possible to treat the 
nuclei as sources of an external field and to calculate the 
probability of n-pair production P n (p) in collisions of two 
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Fig. 1. Feynman diagram for e + e~ pair production in heavy- 
ion collision (first Born approximation). Bold lines denote nu- 
clei, thin lines denote electrons. 
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Fig. 2. Feynman diagram for the unitarity correction to e + e 
pair production in heavy-ion collisions. 



improved approximations for these functions. Thus, the 
aims of the present paper are: (i) to revise the problem 
of the impact-parameter dependent pair production prob- 
ability n e (p), (ii) to update the unitarity corrections for 
the previously calculated processes §5$ and Jl]), (iii) to 
calculate the unitarity correction for photon emission at 
nuclear collisions 



Zi 



Z\ + Z 2 + 7 . 



(9) 



In the latter case, the probability for the production 
of one and only one photon in a heavy-ion collision is 
modified by the necessity of suppressing the possibility 
of simultaneous production of electron-positron pairs. In 
the impact-parameter representation, this implies that the 
cross section for photo-production a 1 should also be mul- 



tiplied by the factor exp( 
replacement 

CT 7 = / P 1 (p)d 2 p — > (7 -y 



-n e ). This corresponds to the 



<r- y +a™ it = J P 7 (p)e-»'W d 2 p, 

where P 7 (p) is the probability to emit a photon in the col- 
lision of two nuclei at a given impact parameter p. There- 
fore, the unitarity correction in this case is given by the 
expression 



nuclei at a fixed impact parameter p [6] . The sum over n 
of the probabilities P n (p) for n-pair production must be 
unity. The unitarity requirement is fulfilled by the Poisson 
distribution 



Pn(p) = ^re- 



(5) 



whose sum over n gives one. Here n e = n e {p) is the aver- 
aged number of produced pairs at a given p, and the factor 
exp(— n e ) is the vacuum-to-vacuum transition probability 



P-t(p) 



1 



,(p) 



d 2 / 



(11) 



It should be mentioned that there are two different 
mechanisms for the photon emission in nuclear collisions: 
the ordinary bremsstrahlung via virtual Compton scatter- 
ing (see Fig. [3]) and the recently considered (see Refs. 
H]) emission of photon via virtual Delbriick scattering as 
illustrated in Fig. [4] 



P (p)=e-^=l-J^P n (p). 



(0) 



tP[ P L m 



,P[ 



Roughly speaking, the probability for producing one pair, 
given in perturbation theory by n e , should be modified 
to read n e exp(— n e ), and this correction is not small for 
an appreciable value of n e . This means that also the cross 
section for the one-pair production <j\ should be multiplied 
by an appropriate factor exp(— n e ) in the integral over 
the impact parameter, which corresponds to the following 
replacement, 



ci = / n e (p) d 2 p 



01 



_umt 



where, finally, 



(7) 



_unit 



n e (p) \l - e -"«H d 2 p 



(8) 



is the unitarity correction to the one pair production cross 
section. 

In Refs. [1] and [S] , rather rough approximations were 
used for the function n e (p). Here, our intention is to use 



(a) 



(b) 



Fig. 3. Ordinary nuclear bremsstrahlung is the emission of a 
photon in a nuclear collision via a virtual Compton subprocess. 



Both of the above mentioned corrections to e + e pair 
production and to photoproduction correspond to the ex- 
change of virtual light-by-light scattering interactions be- 
tween the nuclei. Each block of light-by-light scattering 
brings in an additional factor {Z\u Z^oi) 2 in the amplitude 
of the corresponding process, and therefore, such correc- 
tions can be omitted for the scattering of light ions, for 
muon-nucleus or electron-nucleus scattering. We definitely 
need Z to be large for the correction to be appreciable, and 
consequently this paper is focused on heavy-ion collisions. 

We organize our paper as follows. In Sec. [U we briefly 
recall relevant physical parameters for modern heavy-ion 
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Fig. 4. Emission of a photon in a nuclear collision via virtual 
Delbriick scattering in the lowest order of QED. 



machines and discuss new approximation for the function 
n e (p). Section [3] is devoted to the actual calculation of the 
corrections. Specifically, we consider unitarity corrections 
to e + and p + p~ pair production in Secs. l3.ll and l3.21 re- 
spectively, and we reserve the calculation of unitarity cor- 
rections to the production of photons via virtual Compton 
scattering to Sec. 13.31 and via virtual Delbriick scattering 
to Sec. 13.41 Finally, some conclusions are drawn in Sec. 2J 



2 Toward a Revised Representation of the 
Impact- Parameter Dependent Pair 
Production Probability 

Recently, electromagnetic processes in ultra-relativistic 
nuclear collisions have found strong and partially renewed 
interest in numerous papers (see the reviews [2 , 3 and ref- 
erences therein) . Of topical importance are the RHIC col- 
lider and the future LHC Pb-Pb option. It is therefore 
useful to recall the basic physical parameters of these col- 
liders, namely the charge numbers of nuclei Z\ = Zi = Z 
and their Lorentz factors 71 = 72 = 7. These are given in 
Table [1] which is cited here from Ref . [9] . Many of elec- 
tromagnetic particle-production processes are of imminent 
importance for two reasons: they are either "dangerous," 
e.g. in terms of possible beam losses and background, or 
they are by contrast quite useful for monitoring some ex- 
periments at the RHIC and LHC colliders [TO] . 



Table 1. Nuclear charge numbers Z and relativistic 7 
factors for modern heavy-ion machines 



Collider 


Z 


7 


RHIC, Au-Au 


79 


108 


LHC, Pb-Pb 


82 


3000 



To fix the conventions used, we mention that natural 
units with % = c = 1 and with the fine-structure constant 
a w 1/137 are used throughout the text, and we denote 
the electron and muon masses by m and p, respectively. 

As it was mentioned in the introduction, the impor- 
tance of unitarity corrections has been recognized for the 



first time at the process of electron-positron pair produc- 
tion because a large number of real electron-positron pairs 
are typically produced in ultra-relativistic heavy-ion colli- 
sions. In the lowest QED order (Born approximation) this 
process is described by the Feynman diagram of Fig. [JJ 
Let n e = n e (p) be the expected (average) number of pairs 
to be produced in the collision of two nuclei at a given 
impact parameter p. A closed form of the correspond- 
ing expression was obtained in Refs. [TT1[T2] although a 
complete and consistent interpretation of the expressions 
found was only given later in Ref. [T3]. One issue is that 
n e (p) derived in Refs. [TT]fT2] actually requires a further 
regularization, which was implemented in Refs. |14[I15] . 

As evident from Eqs. (JSJ) and (fl~Tj) . the function n e (p) is 
a very important quantity for the evaluation of unitarity 
corrections. However, the obtained close form for n e (p) 
is, in fact, a nine-fold integral and its calculation is very 
laborious. Therefore, a simpler approximate expression for 
n e {p) is very desirable. 

The properties of n e (p) have been studied in detail in 
Refs. [4lll6j. The corresponding functional form for iden- 
tical heavy nuclei in the region p < 7/771 reads 

n e (p, 7, Z) = (Za) 4 F(x, Z) [L - G(x, Z)\ , (12) 



L = ln( 7 2 ), 



m p . 



Here, 7 3> 1 is the usual relativistic factor (see Table 1). 
The analytical expressions for the functions F(x, Z) and 
G(x, Z) have been obtained in Refs. [41IT5] only for large 
values of the impact parameters, 1 <C x — m p < 7 (see 
the Appendix): 



F(x,Z) = — 



56 ln.x — f(Za) 
9^2 



(13) 



F(x,Z)-G(x,Z) = — 



56 (3/2) lux -/(Za) 



In. 



where 



f{Za) = {Zaf £ 



1 



(14) 



^— ( n\n 2 + (Za) 2 ] 

n—l 1 

is the well-known Bethe-Maximon function [in particular, 
/(79 a) = 0.3129, /(82 a) = 0.3318]. On the other hand, 
for calculations related to unitarity corrections we should 
know the functions F(x, Z) and G{x 1 Z) in the range x ~ 
1. 

In Ref. [16] a closed-form expression for the function 
F(x, Z) at intermediate impact parameters in the form 
of five-fold integral has been given. Tables provided in 
Ref. [16] give a very clear numerical picture of the function 
F(x, Z) for x = 0.01 -j- 100 and several important values 
of Z. For the function G(x,Z), the approximation 



G{x,Z) w 1.5 In (x + 1) + 1.2 



(15) 



(independent of Z) has been given in Ref. [16] as a rough 
indicator of the nonlogarithmic (in 7) term in Eq. (|12[) . 

We improve this approximation by using the results 
from a first Born approximation, and a different parame- 
terization has been employed altogether, namely 



n e {p) = {Zaf [A{x) L - B(x)] 



(16) 
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Table 2. Functions A(x) and B(x) as given Eq. (|16[) . calculated using the approach outlined 
previously in Ref. [IT] 



X 


A(x) 


B(x) 


X 




A(x) 


B{x) 


X 




A{x) 


B{x) 


0.01 


3.537 


8.200 


0.40 


1 


028 


2.910 


4.50 





07273 


0.3355 


0.02 


3.082 


7.434 


0.50 





8879 


2.573 


5.00 





06191 


0.2956 


0.03 


2.805 


6.927 


0.60 





7767 


2.285 


5.50 





05327 


0.2619 


0.04 


2.611 


6.535 


0.70 





6895 


2.085 


6.00 





04635 


0.2339 


0.05 


2.460 


6.199 


0.90 





5557 


1.747 


6.50 





04072 


0.2104 


0.06 


2.331 


5.892 


1.00 





5039 


1.614 


7.00 





03606 


0.1903 


0.07 


2.218 


5.613 


1.50 





3293 


1.142 


7.50 





03217 


0.1730 


0.08 


2.122 


5.366 


2.00 





2323 


0.8581 


8.00 





02887 


0.1579 


0.09 


2.037 


5.139 


2.50 





1725 


0.6716 


8.50 





02604 


0.1446 


0.10 


1.962 


4.949 


3.00 





1331 


0.5431 


9.00 





02360 


0.1329 


0.20 


1.486 


3.900 


3.50 





1066 


0.4560 


9.50 





02149 


0.1225 


0.30 


1.217 


3.345 


4.00 





08717 


0.3879 


10.0 





01965 


0.1132 



(for a list of numerical values, see Table [2]). There is an ob- 
vious connection between the above two sets of functions 
in the limit of low nuclear charge numbers, 



F(x, Z -> 0) = A(x) , G{x,Z^0) 



B(x) 
A(x) 



(17) 



First of all, it is reassuring to verify, based on numerical 
data presented in Refs. [TBlfTT] . that the equality F(x, 0) = 
A(x) is valid in the phenomenologically important interval 
x = 0.01 -r 10 with an accuracy better than 5%. 

One can now take the data for F(x, 0) as given in 
Ref. TBj and use a least-squares method in order to fit 
G(x, 0) to numerical data in Table [21 assuming the func- 
tional form (|16p . A least-squares fit, assuming the depen- 
dence 



G(x,0) = 1.5 In (a; + a) + b , 



(18) 



gives the best estimates a = 1.4 and b = 1.9, where the 
prefactor of the logarithm is fixed by the ratio of the two 
quantities discussed in Eq. (|13| . The approximation thus 
obtained differs from B(x)/ A(x) in the phenomenologi- 
cally most important interval x = 0.02 -j- 5 by less than 
2%. 

Based on the deviation of F(x, Z — 0) from F(x, Z) by 
no more than 25% due to Coulomb corrections for heavy 
nuclei (TSJ, we expect that the same deviation is valid at 
comparison of G(x,Z = 0) and G{x,Z). Taking into ac- 
count that the function G(x, Z) is the subleading term of 
the relative order of 1/L, we can conclude that the ap- 
proximate expression 



n e (p) w {Zaf F(x, Z) [L - 1.5 In (a: + 1.4) - 1.9] 
L = In (7 2 ) , x = m p , 



(19) 



which involves the function F(x, Z) from Ref. [IS] , has an 
accuracy of the order of 5 %. In the calculations reported 
below, we use this very expression. 



3 Calculation of Unitarity Corrections 

3.1 Unitarity corrections for the e+e — pair production 

Unitarity corrections for the process have been consid- 
ered in Refs. [4, 18J. Based on Eq. ([8]), we find the unitarity 
correction is 



—unit 

a e+e- 



1 



-n e (p) 



,(p)d 2 p, 



(20) 



where the integration limits (minimum and maximum im- 
pact parameters) are to be specified below. Since we are 
interested in so-called "silent events" without any "touch- 
ing" of the nuclei, the physically allowed minimal value of 
the impact parameter is 



Pmin — 2i? , 



(21) 



where R is the nuclear radius. A priori, the upper limit 
is Pmax = oo, but due to fast convergence of the inte- 
gral we can use as well /O m ax = 100/m for a quantitative 
estimate of the unitarity correction. Indeed, for large p, 
the asymptotics (|13|) allow an expansion of the exponen- 
tial for the vacuum persistence amplitude in (|10p and in 
analogous expressions used for the unitarity corrections in 
this article, and thus there is a sufficiently large negative 
power of p characterizing all integrands for large p; we 
can thus neglect of excessively large impact parameters in 
the evaluation of all unitarity corrections. As a result, the 
integration region in the variable x can be safely chosen 
as 



xq < x = mp < x n 



100, xq = mp-a 



2mR. (22) 



Here, m is the electron mass. Using the relevant physical 
parameters, we find x$ = 0.0361 for Au, a; = 0.0368 for 
Pb, x = 0.0385 for U and x = 0.0213 for Ca. Finally, 
the relative magnitude of the unitarity corrections for the 
considered e + e~ pair production process ([3]) is 



_unit 



(23) 
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where the known Born cross section reads [19,20 
28 

CTBom = <?o [L 3 - 2.198 L 2 + 3.821 L - 1.632] (24) 



with 



{Z 1 aZ 2 a) 2 

cto = 5 , L = ln(7i7 2 ) . (25) 



For light nuclei with not excessive nuclear charge num- 
ber [(Za) 4 L -C 1], it is possible to calculate in the Born 
approximation the following integrals 



C n = — / F n (x)xdx, n = 2,3,4, (26a) 



D 2 = 2n F 2 (x)G(x)xdx, 
Jxo 



(26b) 



(26c) 



where F(x) = F(x, Z = 0) and G(x) = G(x, Z = 0). They 
can be use to calculate, for light nuclei, 

— the unitarity correction for light nuclei, 



unit 
7 e+e- 



_ 2 i-^2! (C 2 L 2 - D 2 L + E 2 ) , (27a) 



— the total cross section for the production of two e + e 
pairs in collisions of light nuclei, 



02 



{Zaf 



77? " 



(<7 2 L 2 -D 2 L + E 2 ) 



(27b) 



— and the leading logarithmic asymptotics for the total 
cross section a n for n-pair production with 77 > 2, 



_{ZaT „. 

"n — o L'n 



(27c) 



Using the parameterization (|19p , we obtain the follow- 
ing numerical results for the coefficients listed in Eq. (|2"6"|) , 
which enter the formulas listed in Eqs. (|27a[) . (|27b[) . and 
(|27c|) for various unitarity corrections and cross sections 
[in all integrals we used the integration region l]22[) with 
xo = 0.0213, which is the value obtained for a typical 
nucleus of low charge number, namely Ca] 

C 2 = 2.21, C 3 = 0.443, (74 = 0.119, (28) 
D 2 = 15.5, £ 2 = 28.9. 

The results from the previous investigations in Ref. [1] 
concern only those coefficients which can be defined ex- 
clusively in terms of F(x). The previous results read 
C 2 = 1.33, C 3 = 0.264, and C 4 = 0.066. These differ 
from the new results listed in (|28p . because in the previous 
investigation, a less accurate representation of F(x) was 
used, which leads to discrepancies especially when higher 
powers of F(x) enter the integrands as given in (f26|) . 



As an example, using the result (|28|) . we found that 
the cross section for production of two e + e~ pair for the 
Ca-Ca collisions at the LHC collider (7 = 3700) is 



cr 2 = 0.114 barn. 



(29) 



For heavier nuclei, one cannot use the Born approxima- 
tion F(x, Z) w F(x, Z = 0) = F(x) any more. In this 
case, one has to resort to numerical data given in the Ta- 
ble 1 of Ref. [16J for the heavy collision systems Au-Au 
and Pb-Pb, and employ the relativistic factors as given 
in Table [TJ Indeed, unitarity corrections for the process 
Z\ Z 2 — > Z\ Z 2 e + e~ have been considered in Ref. [4] and 
estimated as S e + E - = —4.1% for Au-Au at RHIC and 
S e + e - = —3.3% for Pb-Pb at LHC, where exactly the rel- 
ativistic factors as given in Table [T] have been employed. 
We recall that the ratio 5 e + e - has been defined in Eq. (f23|) . 
Using Eq. (|19p , we are now in the position to give the new 
values 



-5.0% for the RHIC, 
-4.0% for the LHC, 



(30) 



which differ from those obtained in Ref. [4] by about 20%. 



3.2 Unitarity corrections for the fi +/x pair production 

Unitarity corrections for the process (j4|) have been roughly 
estimated in Ref. [5]. Based on the considerations leading 
to Eqs. ([8]) and (fTTj) . we can immediately write down the 
corresponding formula as 



1 -e" 



■Xp) 



MP) d 2 



P, 



(31) 



where Pb{p) is the probability to produce muon pair in 
collisions of two nuclei at a given impact parameter p in 
the Born approximation [the Coulomb corrections to this 
probability, which correspond to multiphoton exchange of 
the produced e + e~ with the nuclei, are parametrically 
suppressed due to the large muon mass and can be ne- 
glected — see Ref. [5] for detail] . 

For a simple calculation, we can use an expression 
for Pb(p) given in the leading logarithmic approximation 
(LLA) in Ref. 0. It reads 



MP) 



28 (ZiaZ 2 a) 

— — 773^ — • 



9tt 



(M 2 



(32) 



Depending on the value of p, the function 'P(p) assumes 
two different asymptotic forms, as shown in Ref. 5 , 



Hp) 



4 In ( — 

\PP 



R 



R 



$(p) = In 2 



pfpR 



P~ p 2 R 



R<^p<-, 
P 
(33a) 



(33b) 
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This expression is valid for large values of ln(p/R), that is 
correct for LHC but not for RHIC. Therefore, we consider 
below the case of the LHC collider only. 
Using formulae (f3T)) -([3"3" ]) we obtain 



O'u+u- 



-49% for the LHC . 



(34) 



where 3^+^- is of course defined as the relative magnitude 
of the unitarity correction in comparison to the Born cross 
section 0^+^-, in analogy with ([25)1 . 

The roughly tenfold increase of the unitarity correc- 
tion ()34jl for muon pair production in comparison to (|30j) 
for electron-positron pair production demands a qualita- 
tive explanation. Indeed, the importance of the unitarity 
correction is due to the enhanced contribution of the re- 
gion of small impact parameters in the impact-parameter 
dependent muon pair production probability (|32[) . Due to 
the prefactor 1 / p 2 in ([32)) , the unitarity correction is log- 
arithmically enhanced as it involves an integration pro- 
portional to J d 2 p/ p 2 over the range of impact parameter 
2R< p < 1/m. 



3.3 Unitarity corrections for Compton-type 
photoproduction 

To a good approximation, tree-level photon emission by 
nuclear bremsstrahlung is described by the block Feynman 
diagrams of Fig. [3J Let the cross section da^ T and der£ r cor- 
respond to the diagrams of Fig. (3^a) and Fig. (3]Jb) , respec- 
tively. Roughly speaking, diagram (a) describes the emis- 
sion of radiation by the Compton scattering of an equiva- 
lent photon, generated by nucleus 1, off nucleus 2, whereas 
for diagram (b), the situation is reversed. The cross sec- 
tion for photoproduction by nuclear bremsstrahlung then 
is obtained as the sum 



der br 



dcr br ! 



(35) 



because the interference term is small and can be safely 
neglected. 

In the LLA, the cross section dcr£ r can be calculated 
using the equivalent photon approximation, in which it is 
expressed as follows: 



derg r = dni dcrc(w, E 7 , E 2 ,Z 2 ) 



(36) 



Here, dni is the number of equivalent photons emitted by 
nucleus 1 in the energy interval duj and the impact param- 
eter range d 2 p, and d<7c(w, -Ely, E 2 , Z 2 ) is the differential 
cross section for the Compton scattering off nucleus 2, for 
an energy E-y of the emitted photon, and an energy E 2 
of the second nucleus of charge number Z 2 and mass M 2 . 
The number of equivalent photons reads 



dni = 



Z\ a dtu d 2 p 



(37) 



and 



E~, 



4 7 f (1 - x 7 ) 



(39) 



For the Compton cross section, we can use the follow- 
ing well-known expression, which is valid for a nucleus 
approximated by a charged point particle. This approach 
gives a good approximation at least in the region of not 
too energetic photons, where the nuclear structure can be 
safely neglected, and it reads 



da c (u,E 2 ,E 7 ,Z 2 ) 



(40) 



4^,2 r 



47T • 



Z 2 ot 



(1 - x 7 ) (y - 2y 2 + 2y 3 ) + 



12 

yx 2 



dRy 

E~, 



where 



Ej 



y 



(41) 



Then we integrate Eq. (|36[) over to and write the result in 
the form 

dag r = dP a (p)d 2 p, (42) 
where the differential probability dP a (p) assumes the form 



dPa(p) 



Z\a a^(Z 2 



1 



4*7 



dE 7 

E~, 



with the Thomson cross section 
ctt{Z 2 ) 



nZla 2 



(43) 



(44) 



Formally the probability in Eq. (|43[) is divergent if inte- 
grated over all photon energies. However, using reasonable 
upper and lower bounds for Ej, the probability remains 
small even at the minimal impact parameter. 

According to the parameter region relevant for the 
equivalent-photon approximation (|37p . this expression is 
valid in the dominant region 2R < p < 7i/w m i n . Integrat- 
ing Eq. (|42p over this region, we obtain 



d< r = 2 



Z 2 a 



a T (Z 2 ) 1 



3 \ dE, 

4 X t 1 E^ 



where 



L = m ( fag.) = l n 
7 V 2R ) 



f 2 7i7l (! ~ x i) 

\ R E-y 



(45) 



(46) 



Now the unitarity correction <5 7 , expressed as a fraction 
of the complete nuclear bremsstrahlung cross section, can 
be obtained by considering diagram (a) alone, 



r <Knit 



(47) 



with the integration region 



< lj < 



7i 



2i? < P < Pn 



7i 



(38) 



and it can be calculated using 



dcr unit = - 



dPa(p) 



1 -e 



-n e (p) 



d 2 / 



(48) 
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The main (logarithmically enhanced) contribution to 
dcr" nit is given by the impact parameter region 2R < p < 
mT 1 , and therefore, a simple estimate can be given as 



-{Zaf— In 



1 



2Rm 



(49) 



A more accurate calculation is based on the direct inte- 
gration of the vacuum persistence amplitude against pair 
production that involves the number of produced electron- 
positron pairs according to Eq. (|12p . and reads 



<5 7 = 



1 



100/m 



J 1 J2R 



1 



-n e (p) 



dp 
P 



(50) 



where the convergence of the integral is assured by the 
asymptotics given in Eq. (|13p and the upper limit of 
100/m for p could have been replaced by oo. An eval- 
uation based on Eq. ([50)) gives the following result for 



P 7 = 1 GeV, 



-19% for the RHIC , 
-15% for the LHC . 



(51) 



3.4 Unitarity correction for Delbriick-type 
photoproduction 

Photon emission in heavy-ion collisions via virtual 
Delbruck scattering has recently been considered in 
Refs. p3|5], where a surprisingly large cross sections was 
found for this case: 



ct 7 d = 14 barn for the RHIC , 
(t 7 d = 50 barn for the LHC . 



(52) 



The main contribution to these cross section comes from 
the photon-energy region 



m <C E.y -C 7m . 



(53) 



Below, we estimate the unitarity correction for the dis- 
cussed process. 

The probability P 7 d (p) , entered the cross section <jTTJ) 
for the photoproduction via the virtual Delbruck scatter- 
ing, can be easily obtained as a function of the impact pa- 
rameter p in the dominant range 1 <C rap <C j 2 , but the 
unitarity correction is mainly given by integration in the 
region mp ~ 1. Therefore, while we estimate the unitarity 
correction here, we stress that a more accurate calculation 
would require a direct evaluation of the Delbriick-type 
photoproduction probability P 7 d(p) in the range mp ~ 1, 
which was beyond our scope. 

For the sake of simplicity, we consider here only the 
symmetric case 7 = 71 = 72 with identical nuclei Z = 
Z\ = Zi. The cross section <t 7 d can be expressed by the 
integration of P 7 d(p) over the impact parameter via the 
relation 

<x 7 d - / P 7D (p)dV (54) 



In LLA we can use the differential cross section in the 
same form as in (|36|): 

(55) 

where <7d(Z) is a high-energy limit of the Delbruck scat- 
tering cross section defined according to Eq. (7) of Ref. [5] 
and 

Z 2 a dui d 2 p 

— r ' 56 

p 2 



der 7 D = 2dn 7 an(Z 2 ) 



is the number of the equivalent photons. The factor 2 in 
right-hand-side of dn 7 takes into account two possibilities, 
corresponding two diagrams in Fig. 2] This expression is 
valid in a parameter range satisfying the two conditions 
m/7 w ^5 ll P an d ^ P ^ 7 2 /m. After integration 
(I56[) over u>, we obtain the probability P 7 d(p) in the form 

2Z 2 a av(Z) ^ ( 

mp, 



PMP) 



(57) 



Under the restrictions for which the approximation made 
in Eq. ([5*6")) remains valid, the expression ([57]) is applicable 
in the dominant region 

2 



1< mp < 7 



(58) 



where again the probability is small. We rewrite (fTTj) in 
the form 



-7D 



pMp) 



1 



d 2 



P, 



(59) 



and now we can use this formula in order to estimate the 
relative magnitude of the unitarity correction. The func- 
tion n e (p) is of the order of unity at p ~ 1/m and given 
by the expression 



n e (p) w 0.5 (Za) 4 L 



(p w 1/m) , (60) 



but drops very quickly at larger impact parameters, with 
an asymptotic behavior of n e (p) oc l/p 2 . Since the func- 
tion P 7 d(p) also drops at large impact parameters [see 
Eq. (|57p ], the main contribution to cr^ 1 ' comes from the 
region p ~ 1/m, and we can estimate the integral (|59[) as 
follows, 



' 7 d 



P 7 d(/o) n e (p) d 2 p 



-P 7 D (1/m) n e (1/m) . 



(61) 



Taking into account the result (|57[) . we obtain the estimate 

at p~l/m. (62) 



/ x 2Z 2 a , N , 
PMp) 5- <td(Z) m 2 L 



The relative magnitude of the unitarity correction thus is 

,-unit 

(63) 



<5 7 d = ^ 0.5 {Zaf . 

<7 7 D 



For the nuclear collisions at modern heavy-ion machines 
with parameters as listed in Table [TJ one can estimate uni- 
tarity corrections to the photon emission to be on the level 
of —5%. Let us emphasize, that the form of this correction 
(|63p is different from that for the lepton pair production 
discussed in Sees. 13.11 and 13.21 
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4 Conclusions 

In this article, we have considered unitarity corrections for 
e + e~ and p + p~ production in heavy-ion collisions, and for 
the production of photons by nuclear bremsstrahlung and 
by virtual Delbriick scattering. The main results of the 
current investigation can be found in Eqs. (|27a|) for the 
unitarity correction to e + e~ production for collisions of 
light nuclei, in Eq. (|30p for the same process in heavy-ion 
collisions in modern colliders with parameters as detailed 
in Table [IJ for /J, + fi" collisions in modern colliders [see 
Eq. (|34[) ], where the unitarity correction is numerically 
large, and in Eqs. (jSTj) and (|63|) for the unitarity correc- 
tion to photoproduction in heavy-ion machines, with al- 
lowance for both the ordinary nuclear bremsstrahlung and 
the virtual Delbriick scattering process. 

Our results as presented for electron-positron pair pro- 
duction in Eqs. (|27a|) and ([30|) are based on a refined treat- 
ment of the vacuum persistence amplitude against multi- 
pair production implied by Eq. (p~9|) . and they represent 
an update of results previously presented in Ref. [4] for 
the same corrections. For production, we update 

the results of Ref. [5]. For ordinary (Compton-type) and 
Delbriick-type photoproduction, the results for the uni- 
tarity corrections are obtained here for the first time to 
the best of our knowledge. 

Finally, we notice that the estimates given here for 
the coefficients C, D and E in Eq. (|28|) also enter the 
total cross section for the production of two e + e~ pairs 
in collisions of light nuclei [see Eq. (|27bj) ] and the lead- 
ing logarithmic asymptotics for the total cross section a n 
for n-pair production with n > 2 [see Eq. (|27cl) ]. From 
a phenomenological point of view, it is important to re- 
mark that all unitarity corrections reduce the one-photon 
or one-pair production cross sections, and that they can 
be numerically large [see Eqs. (|34| and (f5"l"|) ]. 
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Appendix 

In this appendix we briefly recall some details regarding 
the derivation of Eq. (TnT|) . The functions F and G from 
this equation enter the cross section of the process (j3|) as 



follows, 

d<7i = n e {p) d 2 p = (Zia^a) 2 F(x, Z)[L — G(x, Z)\ d 2 p . 

(64) 

Beyond the Born contribution, we should take into ac- 
count, for this cross section, the so-called Coulomb correc- 
tions with a multiphoton exchange between the produced 
pair and the first or second nucleus: 

dai = d(7Born + ddCoul ■ (65) 

The Born contribution to the functions F and G was con- 
sidered in detail in Ref. [2J. Thus, we only need to find 
the Coulomb corrections which enter Eq. (I13|) as items 
proportional to f(Za). 

The Coulomb contribution to the total pair cross sec- 
tion was calculated in |21] ; it can be presented in the form 

dcr C oui = dni cr Co ui(7i^2 -» e + e~Z 2 ) 

+ dn 2 ct C oui(72-Zi -> e + e~Zi) , (66) 

where 

Zfa dut d 2 p 

dm = — ^ T (67) 

tH uji p 1 

is the number of the equivalent photons, produced by the 
ith nucleus, and 

*Coui(lZ ^e+e-Z) = ~^-f(Za) (68) 

is the Coulomb correction to the total cross section of the 
photoproduction 7Z — > e + e~ Z taken from the well-known 
Bethe-Maximon formula. Integrating Eq. (|67[) over u>i in 
the main region 

772 

< w i,2 < 71,2 m, (69) 

72,1 

we find 

dacoui = (70) 
28 (Z 1 aZ 2 a) 2 ( Hl2\ , 2 

97r z ympy \ mp J 

Comparing this expression with Eqs. and (|rj5|) and as- 
suming Z\ = Z 2 = Z, we obtain both terms proportional 
to f(Za) in Eq. JT3J) . 
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